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LONG-TERM  GOALS 

My  long-term  goal  is  a  complete  and  thorough  understanding  of  the  properties  of  acoustic  pulses  sent 
over  basin  to  global  scales.  In  particular,  I  want  to  understand  the  forward  problem  for  calculating 
travel  times  of  the  early  ray  arrivals  in  long-range  acoustic  transmissions  and  to  understand  the 
sampling  associated  with  those  arrivals. 

OBJECTIVES 

This  work  aims  to  determine  the  extent  to  which  existing  models  of  ocean  variability  can  be  used  for  the 
study  of  long-range  acoustics.  To  accomplish  this  goal,  new  tools  were  developed  to  manage  the  often 
large  size  of  the  model  output,  to  extract  and  construct  the  relevant  acoustic  properties  (e.g.,  full-depth 
sections  of  sound  speed)  from  the  model  output,  and  to  make  the  acoustic  calculations.  Another 
objective  was  to  examine  data  obtained  on  deep  hydrophone  arrays  during  the  SPICEX  experiment  to 
establish  general  properties  of  receptions  that  occur  in  the  shadow  zone. 

APPROACH 

This  project  consisted  of  two  separate,  but  not  entirely  unrelated,  investigations.  First,  long-range 
acoustic  data  obtained  in  the  North  Pacific  during  LOAPEX  (Long-range  Ocean  Acoustic  Propagation 
Experiment)  and  SPICEX  (Figure  1)  were  to  be  analyzed  to  develop  a  quantified,  phenomenological 
description  of  stable  “ray  like”  arrivals  measured  by  deep  hydrophone  arrays.  These  arrivals  appear  at 
travel  times  associated  with  the  lower  cusps  of  the  acoustic  time  front  predicted  by  ray  calculations,  but 
the  depth  of  the  receiver  lies  well  below  the  depths  of  the  predicted  cusps  (Dushaw  et  al.  1999)  (Figure 
2).  The  phenomenological  description  of  these  “shadow-zone”  arrivals  was  to  be  used  to  test  the  results 
from  theorists  and  numerical  modelers  striving  to  explain  the  origin,  and  calculate  the  properties,  of 
these  arrivals.  Second,  state  estimates  from  high-resolution  ocean  models  for  the  North  Pacific,  e.g.,  3- 
D  fields  of  temperature  and  salinity,  were  to  be  used  to  examine  their  suitability  for  making  accurate 
long-range  acoustic  calculations. 

Given  the  difficulty  and  expense  of  accurately  characterizing  the  ocean  environment  for  acoustics  by 
data,  model  state  estimates  may  eventually  provide  a  way  to  obtain  time-dependent  acoustic 
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Figure  1.  The  2004  NPAL  Array. 


[The  2004  NPAL  Array  consisted  of  three  SPICEX  moorings  labeled  SI,  S2,  and 
VLA,  augmented  by  the  75-Hz  Kauai  acoustic  source  and  bottom-mounted  receivers 

denoted  f,  k,  n,  o,  p,  and  r.J 


Figure  2.  Shadow  zone  receptions  from  NPAL  acoustic  source  SI  to  bottom-mounted 

receiver  O. 

[Shadow  zone  receptions  from  transmissions  from  NPAL  HLF-5  acoustic  source 
SI  (250  Hz),  located  near  the  sound  channel  axis,  to  bottom-mounted  receiver  O. 
Receptions  over  the  course  of  two  days  are  shown  in  the  waterfall  plot  with  6 
receptions  per  day.  The  prediction  for  the  acoustic  time  front  using  rays  is  shown  in 
the  upper  right,  with  the  receiver  depth  denoted  by  the  red  line.] 


environments  for  acoustic  studies,  or  for  accurate  acoustic  predictions,  over  long  ranges.  Since 
numerical  ocean  models  are  at  present  mesoscale  resolving  at  best,  the  effects  of  internal  waves  will  be 
modeled  separately  and  combined  with  the  ocean  model  state  estimates  to  approximate  the  most 
realistic  acoustic  environment. 

WORK  COMPLETED 

M  work  has  continued  to  focus  on  the  acoustical  calculations  in  numerical  ocean  models;  a  manuscript 
based  on  this  work  was  submitted  for  publication  in  Fall  2008.  Several  models  for  the  temperature  and 
salinity  in  the  North  Pacific  Ocean  have  been  processed  and  developed  so  that  ray  tracing  in  them  is 
routine.  Figure  3  shows  sound  speed  fields  at  300  m  calculated  for  four  estimates  of  the  ocean  state  in 
midwinter  2003.  From  simple  to  sophisticated,  these  models  are: 

The  2005  World  Ocean  Atlas  (Locamini  et  al.  2006,  Antonov  et  al.  2006).  The  world  ocean  atlas  is  an 
important  element  of  this  work,  boring  though  it  may  be,  because  acoustic  predictions  using  this  atlas 
are  known  to  be  a  reasonable  base  state.  Acoustic  arrival  patterns  calculated  with  this  atlas  are 
reasonably  accurate  in  their  absolute  travel  time  and  dispersal  of  the  multipath  travel  times.  The  atlas 
thus  provides  an  important  test  for  the  acoustic  accuracy  of  more  sophisticated  numerical  ocean  models. 
Sound  speeds  derived  from  the  atlas  can  also  be  used  to  correct  the  time-mean  state  of  sound  speed  in 
the  models;  often  the  time-mean  state  of  sound  speed  from  a  model  has  unphysical  properties. 

Figure  4  shows  a  comparison  of  data  obtained  on  the  Kauai  source  to  receiver  k  (see  Figure  1 )  acoustic 
path  with  ray  predictions  using  the  monthly  realizations  of  the  atlas.  For  this  work,  the  ray  code 
employed  is  that  of  Dushaw  and  Colosi  (1998). 

Smoothed  estimates  of  temperature  derived  from  available  hydrography  and  altimetry  (Willis  et  al. 
2003,  Willis  et  al.  2004).  Willis  et  al.  have  calculated  smoothed  estimates  of  upper-ocean  world  ocean 
temperatures  using  hydrography  (e.g.,  Argo  float  data)  and  altimetry  using  objective  mapping 
techniques.  These  estimates  have  coarse  resolution  and  are  smoothed  temporally  with  a  1-year  running 
mean.  Thus,  this  product  does  not  include  internal  wave,  mesoscale,  or  seasonal  variability.  The 
upper-ocean  temperature  profiles  have  been  extended  to  the  deep  ocean  to  allow  acoustic  calculations 
using  cubic  spline  techniques  and  the  world  ocean  atlas  for  the  abyssal  values.  Salinity  estimates  were 
also  not  yet  available,  salinity  values  were  assigned  using  T-S  relations  from  the  World  Ocean  Atlas; 
salinity  values  are  required  to  calculate  sound  speed. 

The  “ Estimating  the  Circulation  and  Climate  of  the  Ocean  ”  (ECCO)  global  ocean  model  (Marshall  et 
al.  1997a,  1997b,  http://www.ecco-group.orgA.  The  ECCO  model  is  a  data  assimilating  model  that 
attempts  to  incorporate  all  available  data  (e.g.,  altimetry,  Argo  float  data)  to  estimate  the  ocean  state. 
There  are  various  flavors  of  ECCO  ocean  models;  I  am  using  the  1 -degree  resolution,  46-layer  model 
from  JPL  that  uses  a  Kalman  filter  for  data  assimilation.  Newer  versions  of  this  model  have  greater 
resolution,  hence  provide  more  realistic  estimates  of  the  ocean  state  for  acoustics  (see  below  “Perth  to 
Bermuda”.  The  1 -degree  model  suffers  from  a  number  of  problems,  one  of  which  is  unphysical  sound 
speed  gradients  in  the  time-mean  state.  These  gradients  make  the  results  of  acoustic  calculations  rather 
unphysical  (if  not  unrecognizable),  so  the  time-mean  state  of  the  model  is  replaced  by  the  World  Ocean 
Atlas  (with  the  model  contributing  estimates  of  the  variability).  I  expect  that  any  existing  model 
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Figure  3.  A  comparision  of  four  sound  speed  fields  of  the  eastern  North  Pacific  derived 

from  different  data  and  analysis. 

[Sound  speed  at  300-m  depth  derived  from  four  realizations  of  the  ocean:  (1)  the  2005 
World  Ocean  Atlas ,  (2)  objective  maps  by  Willis,  (3)  the  ECCO  model,  and  (4)  the  POP 
model.  The  three  bottom  realizations  are  all  obtained  for  midwinter  2003.  The  top 
panels  show  sound  speed,  while  in  the  bottom  panels  the  world  ocean  atlas  sound 
speed  has  been  subtracted  from  the  model  sound  speed  to  show  the  sound  speed 
“anomaly".] 

will  have  one  pathology  or  another  in  terms  of  its  acoustic  properties;  working  around  such  issues  is 
one  of  the  challenges  of  this  project. 

A  high-resolution  “Parallel  Ocean  Program  ”  (POP)  global  ocean  model  (Maltrud  and  McClean  2005). 
J.  McClean  provided  me  with  model  output  from  her  high-resolution  POP  model  runs  for  the  North 
Pacific.  While  this  model  is  not  data  assimilating,  it  is  driven  by  estimates  of  actual  ocean  forcing  such 
as  precipitation  and  winds,  and  it  is  mesoscale  resolving.  The  0.1 -degree  horizontal  resolution  and  40 
depth  levels  of  the  model  are  ideal  for  acoustical  calculation  purposes.  The  model  output  is  available 
monthly  but  continues  only  through  2003. 

The  ray  travel  times  calculated  in  the  four  models  are  compared  to  measured  travel  times  in  Figures  4 
and  5.  The  calculated  rays  are  roughly  equivalent  to  the  rays  identified  in  the  actual  receptions. 
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Figure  4.  Measured  and  Model  Travel  Times  for  the  4-Mm  Acoustic  Path  from  the  Kauai 
source  to  bottom-mounted  receiver  k  in  the  central  North  Pacific. 

[Blue  lines:  the  measured  ray  travel  time  time  series.  Gray  lines:  equivalent  travel 
times  calculated  using  the  indicated  ocean  model.  An  equivalent  temperature  scale  is 
given  at  right.  The  red  line  in  the  top  panel  shows  the  expected  trend  corresponding  to 
5  m°C/yr  warming  near  the  sound  channel  axis.  This  nominal  number  for  oceanic 
climate  change  was  originally  suggested  by  Munk  and  Forbes  (1989).  ] 


Of  all  the  comparisons  that  have  been  made,  including  the  Willis  objective  map  based  on  Argo  and 
altimetry  and  the  ECCO  data  assimilation  model,  the  POP  model  shows  the  most  similarity  to  the 
observations.  It  seems  clear  that  the  acoustic  data  type  offers  unique  information  for  testing  or 
constraining  the  basin-scale  temperature  variations  of  numerical  ocean  models. 


Another  look  at  the  I960  Perth  to  Bermuda  antipodal  acoustic  propagation  test.  In  preparing  a 
review  talk  on  long-range  acoustic  propagation  employed  for  remote  sensing  of  ocean  temperature,  I 
had  occasion  to  review  the  work  published  some  20  years  ago  describing  issues  associated  with  the 
acoustic  propagation  over  antipodal  distances,  based  on  a  test  in  1960  (American  Geophysical  Union 
1960;  Shockley  et  al.  1982,  Munk  et  al.  1988,  Heaney  et  al.  1991;  see  also  the  new  web  page 
documenting  the  1991  Heard  Island  Feasibility  Test: 

http://909ers.anl.washington.edu/~dushaw/heard/index.shtml ).  The  1 960  experiment  was  to  test  the 
1940’s  conjecture  of  Ewing  and  Worzel  (1948)  that  sound  could  easily  travel  great  distances  through 


POP  Model 


1998  1999  2000  2001  2002  2003  2004  2005  2006 


1000  1500  2000  2500 

Yearday  1996 


3000 


3500 


Figure  5.  Measured  and  Model  Travel  Times  for  the  4-Mm  Acoustic  Path  from  the  Kauai 
source  to  bottom-mounted  receiver  f  off  the  coast  of  California. 

[Blue  lines:  the  measured  ray  travel  time  time  series.  Gray  lines:  equivalent 
travel  times  calculated  using  the  indicated  ocean  model.  An  equivalent  temperature 
scale  is  given  at  right.  The  red  line  in  the  top  panel  shows  the  expected  trend 
corresponding  to  5  m°C/yr  warming  near  the  sound  channel  axis.  ] 


the  ocean.  The  test  was  a  success  in  that  the  sound  from  explosive  sources  deployed  off  Perth, 
Australia  was  detected  at  the  antipode  near  Bermuda.  A  closer  look  at  the  propagation  paths,  however, 
found  that  if  the  horizontal  refraction  of  the  acoustic  path  was  taken  into  account,  the  sound  should 
have  been  blocked  by  the  continent  of  Africa  (Munk  et  al.  1988).  Heaney  et  al.  (1991)  appeared 
toresolve  the  issue,  by  suggesting  that  the  sound  reflected  off  of  bathymetric  features  to  arrive  at 
Bermuda,  but,  from  a  modem  perspective,  the  data  bases  of  bathymetry  and  sound  speed  employed  for 
their  calculation  were  too  primitive  for  that  explanation  to  be  definitive  or  convincing. 


1  reattempted  the  antipodal-path  calculation  using  the  2005  World  Ocean  Atlas  to  model  the  global 
sound  speed  field.  I  searched  for  direct  (non-bottom  interacting)  acoustic  paths,  but  successful  paths 
were  not  obtained  (Dushaw  2008).  Detailed  sound  speed  fields  derived  from  one  of  the  more  recent 
global  high-resolution  ocean  models  (ECC02,  Figure  6),  provide  some  interesting  suggestions  that  may 
account  for  the  successful  antipodal  acoustic  test.  In  particular,  the  Agulhas  Rings  to  the  east  of  the 
Cape  of  Good  Hope  act  as  quite  intense  refractive  lenses  for  acoustic  propagation,  with  sound  speed 
deviations  of  about  25  m/s  larger  than  ambient  sound  speed.  These  features  may  well  have 
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Figure  6.  The  geodesic  path  connecting  the  location  of  the  1960  explosive  shots  off 
Perth,  Australia  and  the  location  of  the  receivers  south  of  Bermuda. 

[A  high-resolution  ocean  model  snapshot  of  sound 
d  speed  at  300-m  depth  derived  from  the  JPL-ECC02  numerical  ocean  model  shows 
some  of  the  challenges  of  calculating  acoustic  propagation  over  antipodal  distances.] 

provided  the  bending  required  for  the  1960  acoustic  paths  to  successfully  skirt  the  southern  tip  of 
Africa  on  their  way  to  Bermuda. 

Development  of  a  Parabolic  Equation  Propagation  Code.  One  of  the  tools  required  for  understanding 
long-range  acoustic  propagation,  and  the  nature  of  “shadow-zone  arrivals”  in  particular,  is  the  parabolic 
wave  equation  (PE)  method  for  calculating  acoustic  propagation  for  a  given  ocean  section.  This 
solution  for  the  acoustic  propagation  is  “applicable  in  strongly  range-dependent  environments  in  which 
mode  coupling  is  severe”  (Munk  et  al.  1995).  This  past  year,  I  developed  a  FORTRAN  95  version  of 
the  RAM  Parabolic  Equation  code  (Collins  1993a,  1993b)  which  is  parallelized  for  broadband 
calculations  using  the  Message  Passing  Interface  (MPI)  suite  of  subroutines  (Figure  7).  The 
FORTRAN  95  code  was  developed  from  a  Matlab  version  of  the  RAM  code  developed  by  M. 
Dzieciuch,  which,  in  turn,  was  developed  from  the  original  Collins’  FORTRAN  77  version.  The  basic 
aims  of  this  development  were  to  develop  a  code  that  I  was  familiar  with,  that  employed  a  modem 
computer  language,  and  that  was  parallel.  (Interestingly,  FORTRAN  95  is  similar  to  Matlab  script 
code,  so  the  development  was  fairly  straightforward.)  The  FORTRAN  95  code  runs  some  30%  faster 
than  the  Matlab  version  and  is  easy  to  run  in  parallel  mode  on  a  computer  cluster  of  any  size.  A  website 
was  set  up  to  document  the  code  and  make  it  available: 

http://909ers.apl. washinuton.edu/twiki/bin/view/Main/RamFortranCode  .  This  website  complements 
the  similar  pages  describing  the  Matlab  version  of  the  code: 

http://909ers.apl.washinuton.edu/twiki/bin/view/Main/RamMatlabCode  ,  and  the  website  describing  the 
“eigenray”  geometric  ray  code: 

http://909ers.apl.washinuton.edu/twiki/bin/vicw/Main/EiuenRav  .  In  conjunction  with  this  work,  I 
obtained  and  got  working  the  code  for  simulating  internal  wave  variability  presently  employed  at  the 
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Figure  7.  A  comparison  of  a  time  front  calculated  by  geometric  rays  (black  dots)  and 
by  the  Fortran  95  version  of  RAM  developed  for  this  project 
[The  colors  show  the  intensity  of  the  acoustic  pattern  in  dB  as  indicated  by  the  legend 
at  lower  left.  The  ray  “ intensity ”  is  approximated  by  the  density  of  the  ray  arrivals 
along  the  branches  of  the  time  front.  The  late,  deep  arrivals  at  lower  right  are  an 
artifact  of  the  calculation;  they  are  early  arrivals  “wrapped  around"  from  the  left.  ] 

Scripps  Institution  of  Oceanography  (Van  Uffelen,  personal  communication,  2008).  This  code,  based 
on  a  model  developed  by  J.  Colosi  and  M.  Brown,  simulates  snapshots  of  stochastic  internal  wave 
variability  based  on  the  Garrett-Munk  spectrum  for  internal  waves.  The  next  step  is  to  combine  these 
internal  wave  model  snapshots  with  numerical  ocean  model  state  estimates  to  obtain  realistic  estimates 
of  long-range  acoustic  propagation  using  the  PE  model. 

RESULTS 

Significant  effort  has  been  directed  at  manipulating  the  oftentimes  large  model  data  files,  and  merging 
those  data  with  acoustic  propagation  capabilities.  It  is  clear  that  the  model  results  are  not  yet  directly 
suitable  for  acoustic  calculations.  Constraints  preserving  essential  acoustic  properties  might  be  a 
useful  consideration  in  future  model  development. 

Acoustic  arrivals  in  the  shadow  zone  were  ubiquitous  in  the  data  obtained  from  deep  receivers  during 
the  SPICEX/LOAPEX  experiment.  This  rich  data  set  looks  as  if  it  will  be  very  useful  for  deriving 
general  properties  of  this  type  of  acoustic  arrival. 


The  acoustic  ray  tracing  code  of  Dushaw  and  Colosi  (1998)  has  been  developed  into  a  subroutine  that 
can  be  called  by  any  computer  code,  such  as  a  numerical  ocean  model  during  its  run.  For  tomography 
and  data  assimilation,  this  tool  is  one  of  the  holy  grails  that  has  not  been  available  until  now. 

The  comparison  of  measured  travel  times  and  travel  times  calculated  in  the  various  model  estimates  of 
the  ocean  state  show  both  similarities  and  differences.  Acoustic  travel  times  detect  oceanic  variability 
that  is  not  captured  by  the  numerical  models.  The  acoustic  data  type  is  a  unique  observation  of  large- 
scale  temperature  that  can  be  used  to  correct  the  numerical  ocean  models. 

There  is  still  a  mystery  about  the  nature  of  the  antipodal  acoustic  path  traversed  during  the  1960  Perth- 
to-Bermuda  acoustic  propagation  experiment.  A  definitive  explanation  for  the  success  of  this 
experiment  remains  elusive.  However,  high-resolution  numerical  ocean  models  provide  estimates  for 
the  sound  speed  environment  that  offer  promising,  if  not  obvious,  avenues  for  further  investigation. 
The  environmental  details  of  high-resolution  models  highlight  the  challenges  of  accurate  acoustic 
calculations  over  long  ranges  in  some  regions  of  the  world’s  oceans. 

IMPACT/APPLICATIONS 

Data-assimilating  models  will  some  day  have  enough  resolution  and  data  constraint  that  they  can  be 
used  for  accurate,  real-time  predictions  of  acoustic  properties  over  any  ranges  anywhere  in  the  world. 
This  work  aims  to  assess  present  capabilities,  and  perhaps  develop  criteria  or  constraints  for  better 
behaved  acoustics  that  modelers  can  begin  to  implement. 

The  acoustical  properties  of  shadow-zone  arrivals,  including  the  basic  oceanographic  and  acoustical 
physics  that  give  rise  to  them,  are  not  well  understood,  although  recently  Van  Uffelen  et  al.  (2008) 
found  that  ordinary  internal  wave  scattering  could  account  for  this  phenomena.  These  properties  are 
relevant  to  the  design  of  deep  acoustic  observing  systems. 

The  development  of  ray  tracing  code  as  a  computer  code  subroutine  means  that  numerical  ocean  models 
can  now  implement  the  basic  acoustic  calculations  routinely.  Such  capability  is  an  important  step 
toward  the  routine  use  of  long-range  acoustic  data  for  data  assimilation  using  general  circulation 
models. 

The  conclusion  that  the  acoustic  travel  times  offer  a  unique  constraint  of  large-scale  temperatures  for 
the  ocean  models  suggests  that  a  long-range  acoustic  array  might  be  a  useful  component  of  an  ocean 
observing  system,  a  complementary  measurement  type  to  Argo  drifting  floats  and  satellite  altimetry. 
These  results  were  submitted  for  publication  in  Fall  2008  (Dushaw  et  al.  2008). 

The  propagation  of  sound  over  very  long  ranges  in  the  Southern  Ocean  is  still  as  yet  not  altogether 
understood  -  how  did  sound  propagation  around  the  southern  tip  of  Africa  on  the  Perth-to-Bermuda 
path?  This  question  has  implications  for  such  things  as  the  performance  of  the  Comprehensive  Test 
Ban  Treaty  Organization  array  of  monitoring  hydrophone  arrays  in  the  Indian  Ocean. 


RELATED  PROJECTS 


This  project  is  a  contribution  to  the  North  Pacific  Acoustic  Laboratory  (NPAL)  collaboration, 
comprised  of  researchers  from  the  Applied  Physics  Laboratory,  the  Scripps  Institution  of 
Oceanography,  and  the  Massachusetts  Institute  of  Technology,  among  others,  (http://nnal.ucsd.edu/) 
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